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Understanding and manipulating work fluctuations in microscale and nanoscale systems are 
of both fundamental and practical interest. For example, in considering the Jarzynski equality 
(e~ /sw ) = e~ l3AF , a change in the fluctuations of e -3 " may impact on how fast the statistical aver¬ 
age of e _3lv converges towards the theoretical value e~ l3AF , where W is the work, f3 is the inverse 
temperature, and A F is free energy difference between two equilibrium states. Motivated by our 
previous study aiming at the suppression of work fluctuations, here we obtain a principle of mini¬ 
mal work fluctuations. In brief, adiabatic processes as treated in quantum and classical adiabatic 
theorems yield the minimal fluctuations in e -/3M . In the quantum domain, if a system initially pre¬ 
pared at thermal equilibrium is subject to a work protocol but isolated from a bath during the time 
evolution, then a quantum adiabatic process without energy level crossing (or an assisted adiabatic 
process reaching the same final states as in a conventional adiabatic process) yields the minimal 
fluctuations in e~ l3W , where W is the quantum work defined by two energy measurements in the 
beginning and at the end of the process. In the classical domain where the classical work protocol 
is realizable by an adiabatic process, then the classical adiabatic process also yields the minimal 
fluctuations in e~ l3W . Numerical experiments based on a Landau-Zener process confirm our theory 
in the quantum domain, and our theory in the classical domain explains our previous numerical 
findings regarding the suppression of classical work fluctuations [G. Y. Xiao and J. B. Gong, Phys. 
Rev. E 90, 052132 (2014)]. 


PACS numbers: 05.40.-a,05.30.-d,37.90.+j,05.20.-y 


I. INTRODUCTION 

In small systems involving few degrees of freedom, 
thermal fluctuations and quantum fluctuations in work, 
heat, and other quantities can be comparable to their 
ensemble mean values. It is of fundamental interest 
to understand and control these fluctuations. To that 
end, fluctuation theorems (e.g., the Jarzynski equality 
m\ and the Crooks theorem ED constitute founda¬ 
tional results because they offer rigorous relations be¬ 
tween nonequilibrium statistical fluctuations with equi¬ 
librium properties. Remarkably, the existence of these 
fluctuation theorems does not rule out the possibility to 
further manipulate the fluctuations in work and heat. 
This was clearly shown, for example, in Ref. [7]. As in¬ 
dicated by the Jarzynski equality, i.e., (e~^ w ) = e~ l3AF 
(where W is the work done on a system initially prepared 
at thermal equilibrium with inverse temperature /3, A F 
is the free energy difference between the initial and fi¬ 
nal equilibrium states of the same /?), the mean value of 
e~ l3W determined by (3 and A F, but the fluctuations 
in the work statistics, as manifested in the variance of 
e~P w , can still be suppressed by an additional control 
field [7, [5J . Such suppression in the work fluctuations 
may be beneficial in achieving a better convergence of 
the statistical average of e~ l3W towards the theoretical 
value e~P AF . This hints the possibility in boosting the 
performace of the Jarzynski equality in an actual appli¬ 
cation. 
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There is a second reason for us to focus on the work 
fluctuations. In particular, advances in our knowledge of 
thermodynamics applied to small systems have simulated 
studies of efficient energy devices (classical or quantum) 
at the microscale and nanoscale EHEEj. In these devices 
the work fluctuations are shown to be significant and the 
suppression of the work fluctuations is seen to be related 
to the lieat-to-work conversion efficiency mmm - in - 
deed, our recent study m also showed that the heat-to- 
work efficiency in nanoscale Carnot cycles may not reach 
the Carnot efficiency because, in general, such energy 
devices are operating at nonequilibrium conditions and 
understanding the effect of work fluctuations therein is 
of fundamental interest. A natural question then follows. 
Is there a lower bound of the work fluctuations and if yes, 
how to reach the lower bound? 

Here we restrict ourselves to work protocols in ther¬ 
mally isolated systems, i.e., in the absence of a heat bath. 
These processes can be one important step in the cyclic 
operations of a quantum heat engine, such as in Carnot 
cycles and Otto cycles 0EHH3. The energy change of 
the system in these processes just reflects the work done 
to the system. This paper considers the quantum domain 
first, where the quantum work is defined via two energy 
measurements in the beginning and at the end of a pro¬ 
cess. In particular, given a quantum state initially pre¬ 
pared at thermal equilibrium with inverse temperature [3 
(the same initial state preparation was assumed in deriv¬ 
ing the Jarzynski equality), we shall demonstrate that 
a quantum (mechanical) adiabatic process (as defined 
in the celebrated quantum adiabatic theorem m), ° r 
a corresponding controlled process reaching the same fi- 



2 


nal states as in the conventional adiabatic process, yields 
the minimal fluctuations in an exponential form of the 
quantum work, namely, . The implicit assumption 

for this principle to apply is that the ordering of the en¬ 
ergy levels of the initial and final Hamiltonians does not 
change. Thus, compared with any other unitary evolu¬ 
tion to realize a work protocol in a thermally isolated 
system, a quantum (mechanical) adiabatic process, exe¬ 
cuted slowly or executed fast via a control field, achieves 
a lower bound in the fluctuations of e~^ w . We term this 
original finding as the principle of minimal work fluctua¬ 
tions. Numerical experiments based on a Landau-Zener 
model dang confirm our theory. 

The principle of minimal quantum work fluctuations is 
also extended to the classical domain. The motivation of 
this extension is two-fold. First, the proof of the principle 
in the quantum domain clearly indicates the existence of 
a semiclassical analog. In particular, the quantum energy 
level index n, which is important in our quantum proof, 
reminds us of the classical action variable / because in 
semiclassical quantization, n and I are simply related. 
Second, in our previous study aiming at the suppression 
of work fluctuations in a classical system, we observed 
that an optimal control field can dramatically suppress 
the fluctuations in but can never outperform that 

achieved from a classical (mechanical) adiabatic process 
in all the cases studied (within statistical error) [7], In¬ 
deed, we are able to explicitly show, at least for one¬ 
dimensional classical systems, that a classical adiabatic 
process yields the minimal fluctuations in . 

This paper is arranged as follows. In Sec. II, we review 
the definition of quantum work with necessary details. 
In Sec. Ill, we calculate the fluctuations in and 

give a detailed proof of the principle of minimal fluctua¬ 
tions in quantum work. This theory is tested in Sec. IV 
in a finite-time Landau-Zener process, by comparing the 
fluctuations attained using an optimal control approach 
with those obtained from assisted adiabatic passage. In 
Sec. V, we extend our minimal fluctuation principle to 
the classical domain. 


II. QUANTUM WORK AND RELATED 
QUANTITIES BASED ON TWO-TIME ENERGY 
MEASUREMENTS 


Let us consider a quantum system with a Hilbert space 
of dimension N {N can be infinity). The system is sub¬ 
ject to a unitary process described by a time-dependent 
Hamiltonian H(t) starting from t = 0 to t = r. For 
i = 1,..., N, let \ipi) and | ip'f) denote the eigenstates of 
H 0 = H( 0) and H T = H(r ), respectively, with eigenval¬ 
ues Ei and E[. Let U denote the unitary time evolution 
operator associated with the whole process from t = 0 to 
t = T. 

At t = 0, the system is assumed to be prepared in the 
Gibbs state with inverse temperature /3, with its equilib¬ 


rium density matrix given by po- 

N 

po = (i) 

i= 1 


with 


o-0 Ei 


Pi = 




s K e 


-0Ei 


( 2 ) 


where Zq is the quantum partition function. 

At the end of the work protocol, one asks how much 
work has been done to the system. This can be answered 
by considering two energy measurements at t = 0 and 
Indeed, for an arbitrary function f(E, E’) 
involving two energies E and E' associated with Hq and 
H t , one needs to perform such kind of two-time mea¬ 
surements, with the average value of f(E,E r ) given by 


N 


f(E,E')=Y J Pi\ a ij \ 2 f(Ei,Ej), 

i,j =1 


( 3 ) 


where 


aij = (), (4) 

with 

N N 

Ek-i 2 = Ei^i 2 = i- ( 5 ) 

i=i j=i 


Experimentally, pi\aij\ 2 can be interpreted as the prob¬ 
ability that one finds the system in the i-th eigenstate 
of Ho at t = 0 and then in the j-th eigenstate of H T at 

t = T. 

Following this definition based on two-time energy 
measurements, the average work done on an isolated sys¬ 
tem is given by 

n 

(W) = ~E^H = V p, a,, 2 {E' - Ei) 

, ; ( 6 ) 

= Trace[/7/?oE^.ff r ] — Trace[po^o]j 


where po is the density matrix right after the first energy 
measurement at t = 0. For our case or any other ini¬ 
tial state as a completely mixed state in terms of energy 
eigenstates, po = Po- 

We note in passing that with the above-mentioned ini¬ 
tial equilibrium state, another definition of the quantum 
work, i.e., d(W) = Trac e[dHp\ (20^ yields the same ex¬ 
pression. That is, 


tyV) = f Trac e[dH(t)p(t)} dt 
Jo 

= I d{Trace[H (f) p(t)]} 

Jo 

= Trace[iL r p r ] — Trace[HoPo], 


( 7 ) 


where p(t) is the time evolving density matrix and we 
have used Trace[dp(t)IL(f)] = 0 for a thermally isolated 
system under unitary evolution. 
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III. PRINCIPLE OF MINIMAL WORK 
FLUCTUATIONS: QUANTUM SYSTEMS 

In the literature, there is the minimal work princi¬ 
ple m that has been attracting considerable atten¬ 
tion recently. That is, for an isolated system initially 
prepared as a Gibbs equilibrium distribution described 
above [Eq. 0 ], a quantum (mechanical) adiabatic pro¬ 
cess, if implementable, yields a minimized average work 
[see Eq. ([6])] for the system among all possible thermally 
isolated processes that start from Hq and end with H T . 
Of particular interest to studies of nanoscale energy de¬ 
vices, this minimal work principle implies that, in or¬ 
der to minimize the so-called disspated work and max¬ 
imize the heat-to-work conversion efficiency of a quan¬ 
tum heat engine, a thermally isolated step should imple¬ 
ment bare adiabatic processes or accelerated adiabatic 
processes [14]. Both as a fundamental question and a 
practical issue (relevant to understanding the stability 
or reliability of the work output of nanoscale energy de¬ 
vices), here we hope to understand the lower bound in the 
work fluctuations. A digestion of the key element in the 
proof of the minimal work principle, together with the 
Jarzynski equality, stimulates us to consider fluctuations 
in e _ ' SH , rather than in W directly. 

For convenience, when referring to a general adiabatic 
process, we mean both a conventional adiabatic pro¬ 
cess or its alternative version based on engineered time- 
dependent fields, such as those processes named “short¬ 
cuts to adiabaticity” (STA) [22l - f27l . We order the initial 
eigenenergy values as 

Ei < E/2 < • • • < Ejf. (8) 


E[ (eigenvalues of H T ) is the same as £), namely, 

E[ < E' 2 < ■ ■ ■ < E' n . (10) 

This leads to 

e~ 2/3 E'j+i < e~ 2 P E 'i, j = 1,2, • • • , N — 1. (11) 

which will be useful below. 

We are now ready to prove that a quantum adiabatic 
process produces a minimized variance in e ~^ w , as com¬ 
pared with all other thermally isolated processes that 
start from Hq and end with H T . To distinguish from 
work in other cases, we use W to denote the work in a 
quantum adiabatic process and use dj 7 to represent the 
associated transition probabilities from state Ei to state 
E'j [see Eq. 0]. 

The fluctuations of , characterized by the vari¬ 
ance of e~P w , can be quantified by 

(<=-'*) = (( 

= (e~^ w ) - (e~e w ) 2 . 

According to the Jarzynski equality (e~ /3W ) — e~ l3AF , 
the second term of the right-hand side of Eq. © is 
a constant for fixed Hq and H T . Thus, essentially we 
are left with treating the first term only. According to 
Eq. 0, the first term can be evaluated by 

N 

(e -W) = £ (13) 

i,j = l 


e -pw_ {e -p W) f\j 


no'* 


It follows that the thermal excitation probabilities have 
the following ordering 

Pi > P2 > ■ ■ ■ > Pn • (9) 

Assuming that Hq and H T can be connected by a quan¬ 
tum adiabatic process, then we require the ordering of 


Using the identity (summation by parts) 

N N N -1 j 

E bjCj = b N Cj - E ( 6 J+i _ b i) E Cfc ’ ( 14 ) 

2=1 2 =1 2=1 fe = 1 

we obtain from Eq. (131 and Eq. 0 the following: 


N 


0 —2/3 W' 


= E 


Pie 


2PEi 


N 

E 

2=1 


| 2 - 


N-l j 

-2/3£' + 1 _ p-2/9 E'j 


2=1 


; £ 


&ik 


N 


N-l 


= e~ 2 P E 


' n E^ e 


2f3Ei 

= 1 2=1 


E(< 


r 2 P E J+i - p-^ e j 


N 


k =1 
3 




|2 

Uik | • 


k=l 


(15) 


In the same fashion, one can look into the variance in e , which is determined by the ensemble average (e 2 i 3W ') J 
with 


N 


- 2m ) = £« 


2/3 Ei 


i—1 


= e _2/3E 


N 


E e 

2=1 


-2/3 E' n 


an - 


N-l 

E(< 

2=1 


-WEj+i _ p—2/313) 


; )E 


ttik 


N 


' N Y.Pi e ^ Ei 


i= 1 


N-l N 

J- (e- 2 ? E N - e ~ 2 ^) ^2 


k =1 
3 


Pie^2 \ a ik\ 


3 = 1 


2/3 Ei |2 

i—1 k—1 


(16) 
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The difference in the fluctuations of e~^ n between a 
general nonequilibrium process and an adiabatic process 
is given by 

__ JV —1 

(e~ 2f)w ) - (e~ 2f>w ) = - E (e" 2/3 ^ +1 - e 

3 = 1 

(17) 

where 

N j 

= E« e2/3Bl EtM 2 - i“ifci 2 i- ( 18 ) 

2=1 fc=l 


Note now that in a quantum adiabatic process, the state 
populations do not change, with 


\aik\ 2 = Sik- (19) 

Further using Eqs. 0 and 0, we obtain from Eq. ( fl8] ) 


N 


= ^E e/3£! E [m 2 ~ i^i 5 


o 


=i fc=i 


1 J 


,PEi 


> 


Z q o 

e P E j 


E m 2 - 1 


= 1 lk=l 

3 


N 


* 1 ^ 5 : 


Gift 


2=7 + 1 fc = l 


E 

2=1 


E m 2 - 1 


7c=l 


p /3Ej N 3 

6 EEi 

0 i—j +1 fc=l 


Z 9 


^2fc 


H + i] 


= o. 


( 20 ) 


That is, the function S , d ehned above cannot be nega¬ 
tive. Returning to Eq. (|17|, we finally arrive at 


(e~ 2flw ) - ( e~ 2l3w ) > 0 ( 21 ) 

or equivalently, 

a 2 ( e -2 (IW) > ct 2 ( e -2 PW'J (22) 


That is, the fluctuations in e ~ ,3W , as characterized by the 
square variance of er !sw , become minimal if the process 
from Hq to H r is a quantum adiabatic process. 

Together with the previously established minimal work 
principle, one can now conclude that given a thermally 
isolated system prepared initially at equilibrium and later 
driven by Hamiltonian H(t) with specified H 0 at the start 
and H t in the end, then a quantum adiabatic process, if 
implementable, not only yields the minimal average work, 
but also the minimal fluctuations in This is one of 

the main results of this work. Before closing this section, 
a few remarks are in order. 

First, it is necessary to highlight the counter-intuitive 
nature of our minimal work fluctuation principle. Given 
that a quantum adiabatic process minimizes the mean 


work, one may naively think that the corresponding 
ensemble-averaged value of e ~ 2/3W , a key quantity in 
obtaining the variance of e~ ,3W , seems to be maxi¬ 
mized. But our proof indicates the opposite, namely, the 
ensemble-averaged value of e~ 2 ^ w is actually minimized 
in a quantum adiabatic process, and so is the variance of 

e-W. 

Second, throughout the proof, we have assumed that 
the ordering of E[ is the same as Ei to ensure the exis¬ 
tence of an adiabatic process to connect Hq to H T . If this 
prerequisite is not satisfied, then it is unclear what kind 
of process can reach a lower bound of work fluctuations. 
Indeed, there a very slow process not necessarily yields 
very small work fluctuations. By itself this constitutes a 
fascinating topic for future study. 

Third, we have emphasized above that a quantum adi¬ 
abatic process may be executed via a control field by tak¬ 
ing advantage of STA [Zl EU22H24!. Consider, for exam¬ 
ple, one type of STA, which is also called counter-diabatic 
driving )22j or transitionless quantum driving [23| . There 
an additional control Hamiltonian is introduced to drive 
a system. In this case, the work to the system is done 
by two sources: the original protocol to tune the system 
from H 0 to H Tl as well as the additional control field im¬ 
plementing an accelerated adiabatic process. In practice, 
one may wish to set the additional control field achieving 
STA to zero at initial and final times, so that the energy 
measurements of the total Hamiltonian in the beginning 
and at the end still reflect the inherent energy eigenval¬ 
ues of the bare system alone. With that understanding 
in mind, it is now possible to achieve minimal quantum 
work fluctuations using very fast processes! 


IV. PRINCIPLE OF MINIMAL WORK 
FLUCTUATIONS MANIFESTED IN A 
LANDAU-ZENER PROCESS 


Here we use the Landau-Zener model mm to illus¬ 
trate and check the principle of minimal work fluctua¬ 
tions. On the one hand, we use STA to realize a fast adi¬ 
abatic passage in this model and then examine the work 
fluctuations. On the other hand, we use an optimal con¬ 
trol theory (OCT) pB] - f5T] to minimize the work fluctua¬ 
tions, with the variance in e~^ w to be minimized. Some 
necessary details regarding an ensemble-based OCT ap¬ 
proach are presented in Appendix A for completeness. 
According to our principle of minimal work fluctuations, 
the results from our OCT can never beat that from an 
adiabatic process. This fact will be numerically checked 
below. Furthermore, it is curious to see how the results 
from OCT may approach that based on STA. 

The Landau-Zener model Hamiltonian, all in dimen¬ 
sionless units, is assumed to be 

H(t) = Z{t)a z +X 0 a x , (23) 

where a z and cr x are the Pauli matrices. For a time- 










5 


TABLE I. Performance of work fluctuation suppression in a Landau-Zener process, in the absence or presence of control fields 
needed for realizing STA or OCT, characterized by the variance in e ~^ w , using 10 6 two-time energy measurement results. The 
model Hamiltonian is described in (Eq. 231, with Z(t) given in Eq. (251, Z(0) = 1.0, Z(t ) = 3.0, r = 0.0001 (duration of the 
process), Xo = 2.0, and the inverse temperature /3 = 0.1. The numerically found transition probabilities are also presented, 
where |1) and |l') are the ground states of the initial and final Hamiltonian, \2') is the excited state of the final Hamiltonian, 
and U denotes time evolution operator. The obtained values of (e _/3lv ) all agree with the theoretical value e -/3AF « 1.040 
obtained from the Jarzynski equality. 


Process 

(e- pw ) 

a(e~ 0w ) 

!<i'|i/|i>l 2 

| <2 / 1 L/| 1) | 2 

bare system 

1.040 

0.202 

0.9341 

0.0659 

STA 

1.040 

0.134 

1.000 

0 

Optimal control 

1.040 

0.134 

1.000 

0 



FIG. 1. (color online) Time dependence of the amplitude 
of a numerically found OCT field, as compared with that of 
a control field realizing STA in a two-level system (Eq. 231, 
with Z(t) given in Eq. (251, Z( 0) = 1.0, Z(t ) = 3.0, t = 
0.0001 (duration of the process), A'o = 2.0, and the inverse 
temperature set to be /3 = 0.1. All the plotted quantities here 
are scaled and hence in dimensionless units. 


dependent Z(t) and a time-independent Xo, the control 
field achieving STA is given by [25H25] 

H STA (t ) = h(Q 0 /2)<Ty, (24) 

where 0o = arccos (Z(t)/Ro) and Rq = \/Xq + Z(t) 2 . 
As an example, we choose [27] 

Z(t) = Z(0) + [(Z(T)-Z(0)](-) 3 

t t ^ 

15 [Z(r) - Z(0)](-) 4 + 6 [Z(r) Z( 0)](-) 5 

T T 

such that the control field H STA (t) is zero at t = 0 and 
at t = t. In our numerical calculations we choose Z(0 ) = 
1.0 and Z(t ) = 3.0 in dimensionless units. We further 
set r to be as small as 0.0001 (as compared with 
such that the process will be highly non-adiabatic were 
there no control field. 


Parallel to this, the control filed in OCT is assumed to 
be A(t)a y with A(0) = 0 and A(t) = 0. This boundary 
condition for A(t) can be satisfied by introducing an ap¬ 
propriate time profile of a cost function (for details see 
Appendix). It can be checked that the system here is 
fully controllable by a y . That is, any unitary evolution 
operator can be yielded in a finite time by considering a 
control field of the a y type J32]. 

To quantitatively characterize the performance in sup¬ 
pressing the fluctuations in e~^ w in STA and in OCT, we 
randomly sample initial energy eigenstates \^i) accord¬ 
ing to the initial thermal probability distribution (the 
first measurement), then evolve them under the total 
Hamiltonian H = Hq + H STA . Next, we again ran¬ 
domly sample energy eigenstates \ipj) according to its 
probability projected on the final state (the second mea¬ 
surement). Individual values of W are obtained from 
Wi = Ej — E il and the fluctuations in e~^ w are cal¬ 
culated from a(e-f> w ) = \J^ J2iLi( e ~^ Wi ~ { e ~^ w )) 2 i 
where the total number of “trajectories” is chosen to be 
M = 10 6 . 

The variances in e~^ w obtained in the bare system, 
under STA and under OCT are presented in Table |T[ 
along with the ensemble-average of e ~^ w , as well as the 
transition probabilities between initial and final energy 
eigenstates. First of all, it is seen that all the three cases 
yield exactly the same average of e ~^ w , consistent with 
the Jarzynski equality. Secondly, both processes under 
STA and under OCT have suppressed the variance of 
e —f)W. remarkably, the performance of OCT in sup¬ 

pressing the variance of e~ l3W does not beat that of STA 
(producing the same variance here), thus confirming our 
expectation that an adiabatic process yields the lower 
bound for fluctuations in e~^ n . An investigation of the 
transition probabilities also offers more insights. One 
observes that the transition probabilities obtained under 
OCT is actually the same as that obtained under STA. 
That is, an OCT field aiming at minimizing work fluctu¬ 
ations and successfully reaching the lower bound of work 
fluctuations (as characterized by the variance of e~ 3W ) 
tends to reproduce, at t = r, the initial populations on 
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each energy eigenstate. This confirms our theory from 
another angle. 

Lastly, in Fig. 1 we compare the time dependence of 
the control field numerically found in OCT with that 
of the STA control field. It is seen that A(t) found in 
OCT is not the same as in the case of STA. The peak 
field amplitude in the OCT case is smaller. This differ¬ 
ence clearly indicates that there are different solutions 
when it comes to minimize the work fluctuations. In 
our proof above, the crucial requirement to reach the 
lower bound of work fluctuations is a^- = 5j :] . i.e., pop¬ 
ulations on states E[ stay the same as the populations 
on states Ep. what happens during the process does not 
matter. An adiabatic process, which keeps the popula¬ 
tions throughout the process, then offers a lower bound 
in work fluctuations. However, in other processes with 
nonadiabatic transitions during the time evolution but 
still the same final-state populations as in an adiabatic 
process, the lower bound of work fluctuations can still be 
reached. 


V. PRINCIPLE OF MINIMAL WORK 
FLUCTUATIONS: CLASSICAL SYSTEMS 


In our previous study [7j treating classical systems, a 
classical version of STA to realize classical shortcuts to 
adiabaticity 0 and a classical OCT were considered to 
suppression classical work fluctuations. It was found that 
the performance of OCT, as quantified by the suppres¬ 
sion of the variance in e ~ l3W , can at most reach that ob¬ 
tained from classical STA (within statistical error). This 
hints the existence of a lower bound of classical work 
fluctuations, which might have been reached by classi¬ 
cal adiabatic processes. Following the technique we used 
in proving the principle of minimal work fluctuations in 
the quantum domain, here we shall prove an analogous 
principle in the classical domain. 

Consider a time-dependent classical system with 
Hamiltonian H c (q,p, A(f)), where (q,p) represents phase 
space coordinates, and A(t) is a time-dependent parame¬ 
ter evolving from t = 0 to t = r, with Hfi = H c (p , q, A(0)) 
and = H c (p, q, A (r)). In the absence of a bath, the in¬ 
clusive work (33] as a function of initial conditions (po, qo ) 
is given by 

W T (q 0 ,po) = H^[q T (q 0 ,p 0 ),p T (q 0 ,p 0 ),T] - H£(q o ,p o ,0), 

(26) 

where [q T (qo,Po),p T (qoiPo)] are simply the time-evolving 
phase space coordinates staring from (qo,po). Then the 
thermal ensemble average (e~ 2/3w ) needed for calculating 
the variance of e~^ w is given by 


~ 20W ) =JfJJ dqodpodqdp e- 2 ^(«>P)" ff o(9o,Po)] 


( 27 ) 


where po(qo,Po) is the initial Gibbs distribution with in¬ 
verse temperature /?: 

Po(qo,Po) = ^e-^o.Po), (28) 

with 

Z$ = I ^° WO) ^0. (29) 


Next, to have a classical adiabatic process connect¬ 
ing Hq to H^, we further assume that H c (q,p, A(f)) can 
be written as H£(I t ), where (/ t ,0 t ) are the action and 
angle variables associated with H c (q,p, X(t)). That is, 
during a work protocol, the Hamiltonian can be always 
expressed as a function of its instantaneous action vari¬ 
ables only. In particular, at t = 0, Hq = Hq(I 0 ); at t = t. 
H= H^(I T ), where / 0 and I T are the action variables in 
the beginning and at the end. According to the classical 
adiabatic theorem, if the parameter A (t) changes slowly 
as compared with the internal frequency of the system 
(ujt = dH gj ^ ), then the action variable It remains in¬ 
variant during the work protocol. Note that, as our proof 
below shows, what is truly essential is the invariance of 
the action variables at the two boundary times and an 
extra condition (analogous to the no-level-crossing con¬ 
dition in the quantum case) about Hq(I 0 ) and H°(I T ). 
That is, what actually happens for 0 < t < r is not really 
essential in reaching the lower bound of the work fluctua¬ 
tions. Note also that we have assumed that the system’s 
Hamiltonian has only one degree of freedom. Extend¬ 
ing the proof below to systems with more than one de¬ 
grees of freedom is possible, however the prerequisite for 
the principle to be applicable seems to be demanding in 
high-dimensional classical systems. 

Equation (27) can assume a different form if we in¬ 
troduce a canonical transformation from (q,p) to action- 
angle variables (1,9) (only as integration variables), 


(e~^ w )=J J J Jdl 0 d6 0 dld9 e-WWW-WrtpoVo) 

6[I T (I o ,e o )-I}6(6 T (I o ,0o)-e}, 

(30) 

where (/ T , 9 T ) is expressed as a function of the initial con¬ 
ditions (Iq,9 0 ). The explicit dependence of ( I T ,0 T ) upon 
(J O ,0 O ) nright not be spelled out below for convenience. 
The initial thermal density now becomes 

Po(Io) = (31) 

z o 

For later use we also note 


J J dId95[I T (Io, 9 q) — I] S[9 t (Iq, $ o ) ~ 9] 

= J I dI 0 d9 0 S[I T (I 0l Oo) - I} S[9r(Io, Oo) - 0} (32) 

= 1. 


Po(qo,Po) S(q ~ q T ) S(p - Pr) 
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In the quantum case we benefitted from the assump¬ 
tion that the ordering of the final energy eigenvalues E[ 
is the same as that of initial energy eigenvalues E^. This 
motivates us to reorder certain integrals in terms of the 
value of Hq. Specifically, we define a nonconventional 
type of integration as follows: 


the variable / according to H 0 (I ), such that the sum¬ 
mation of f[I(n)]AI over n is executed progressively ac¬ 
cording to the value of H 0 (I). This makes it clear that 

J [dI] HS f{I) = J dl /(/), (35) 


mwiss ( 33 ) 

n= 1 

where I is the integration variable and I(n + 1) — I(n) = 
A I represents an infinitesimal interval in the / space, and 

H^[I(n)}<H^[I(n + 1)], n= 1,2 - - ■ . (34) 

The physical meaning of such a resummation is to order 


namely, the ordering in the summation does not affect 
the final sum. Because of the reordering, we also have 

Po[I(n)\ > Po[I(n + 1)], n = 1,2, ■ ■ • . (36) 


which means that the initial thermal excitation probabil¬ 
ity for I(n) is no less than that for I(n + 1). This fact 
will be useful for our proof below. 

Exploiting the newly defined integral above, Eq. (30) 
can be rewritten as 


(e~^ w ) 


jR jR 

J [dI 0 \H§ e 2/3ff Vo(/o) J d0 o J [dI] HS e ~ 2f>H ZW J dO S[I T - I\6[0 r - 8} 


(37) 


In the above equation the lower and upper limits of the integrals over 0 O or 6 are always 0 and 2n. For the reordered 
integral over [d/o]jj= an d [dl]n^, the lower and upper limits of the integrals are denoted I L and I R : namely, the values 
of Iq that gives rise to the lower limit and upper limit of Hq(I 0 ). Either I L or I R can be +oo. Performing integration 
by parts, we obtain from Eq. (37) 


- 2pW) = f [dI 0 \ H§ e 2 TO ) po(Jo) J d6o j de s[6t _ e] 


(38) 


0 -2/3^(/ h ) _ 


II L 


,-2 PK(I) 


J H o Ji l 


[dl']m S[I T -1'} 


Note that this equation is in parallel with our early quan¬ 
tum result in Eq. (15). 

In the same manner, the ensemble average (e~ 2 ^ w ) 


associated with an adiabatic process is found to be 


( e -2/3W) = j [dIo]H§ e 2pH S (I 0 ) po( j o) j ^ j de5[ 0 T _ 9] 


(39) 




II L 


a -2P H*{I) 


H, 


5 JI L 


Whs 5[/o - /'] 


where W, I T = Iq, and d T represent the values of work and action-angle variables at t = r in an adiabatic process. 
Needless to say, this result is analogous to Eq. (16) in the quantum case. 


The difference in the square variance of e M between a general work protocol and an adiabatic process is then 
given by 


(e- 2 ? w ) - (e- 2 ? w ) 


[ d[-e 2 PK{i) ]h- A (I), 
Ji L 


(40) 













where 


jil / j 

A(/) =J [dloUs e 2/3H ° (Io) po(I 0 ) J d0 o j J dd 5(6 r - 9) j [dI'] H§ S(I T - I') 


(41) 


- J dd S(d T - 6) J [dI'] HS 6{I 0 -/')>• 


Equations (40) and (41) are classical analogs of Eqs. © 
and (fl8|. With the assumption that Hq and can 


be connected by an adiabatic process, the instantaneous 
frequency uit = dH gj^ cannot be zero during the entire 
work protocol with 0 < t < r. Therefore, if for I(n) = 
nAI + I L , Hq[I(u+1)] > Hq[I(ti)], then H^\(I(n + 1)] > 
H^[(I(n)], namely, the ordering of the energy in terms of 
the value of the action variable is the same for Hq and 


r 


H Imagine this is not the case, then at a particular 
time 0 < t < r, we must have H£\(I(n + 1)] = H^[(I(n)} 
and hence u t = 0. With these insights, it is now clear 
that the factor d [—in Eq. (401 cannot be 
negative. 


One can split the integration in Eq. (41) over [dIo\H§ 
into two intervals: from I L to I and from I to I R . This 
leads us to 


A (/) = / [dI 0 ]H S 


e PHZ(I 0 ) 


' II 


zz 


d6 0 \ - dd S(d T -d)+ dd 6[d T -9) [dI'} H§ 6(I r - I') 


/ i L 


[dh 


Oj HS 


o0HZ(I o ) 

Zn 


J ddo J dd 5{d T - d) J [dI'] H§ S(I T - I') 


(42) 


Further noticing that (i) the expression inside the {•} in Eq. (42) is not positive, (ii) the integrals over the S functions 
in the second line of Eq. (42) is not negative, we immediately have the following 


A (I) > f iL [dIo]s S j dd 0 l - j dd S(d r - d) + Jdd S(d r -d) j [<U'] HS S(I t - I') 


ri 
ii L ^ 


e P HS(I) [I* 

J i 

e PH‘(I) rl 


[dlohs J dd Q J dd 6(d r - e) Jjdl'Uc 5{I t — I') 


z, 


0 j I L 


[dIo\ H§ / ddo / dd 6(0 T - 9) 


(43) 


e iw-(i) r i R r r r 1 

[dI 0 \ HS J ddo J dd 6(d T -0) j L [dl'Us S(Ir - I') 

f e PH%(I) r pi 

[dIo\H S J dd 0 + —- J dd Jjdl'Ug 


Z 0 Jil 

e PHg(I) Pi 


Z\ 


0 j I L 


= 0. 


In reaching the last two steps, we have used the fact 

that jj L [dIo\Hjj f ddo d(d T — d)S(I T — I') = 1. Hence, we 
have proved that for a classical adiabatic process with 
the above-assumed conditions (analogous to the no-level- 
crossing assumption in the quantum case), 

<j 2 (e~P w ) > a 2 (e~P W ), (44) 

the principle of minimal work fluctuations in the classical 
domain. 

Before ending this section, we remark that our tech¬ 


nique here can be adapted to prove a classical version 
of the minimal (quantum) work principle obtained in 
Ref. |21| . Considering the importance of the minimal 
work principle and the relevance of classical statistics in 
nanoscale systems with relatively high temperature, in 
Appendix B we indeed give such a proof, which is similar 
to our ideas here proving the principle of minimal work 
fluctuations. 
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VI. CONCLUSION 

We have obtained a general principle regarding the 
minimal work fluctuations for thermally isolated systems 
initially prepared at equilibrium and then subject to a 
work protocol. Specifically, if the initial and final states 
can be connected by an adiabatic process as stated in the 
seminal adiabatic theorem, then the variance of e~P w 
reaches the lower bound, as compared with all other pro¬ 
cesses operating between the same initial and final sys¬ 
tem Hamiltonians. This is true in both the quantum 
and classical domains. It is now clear that an adiabatic 
process yields not only the minimal average work, but 
also the minimal fluctuations in e~^ w . This main result 
represents a somewhat counter-intuitive but fundamen¬ 
tal understanding of work fluctuations. The actual proof 
of our principle of minimal work fluctuations also indi¬ 
cates that the lower bound of work fluctuations can be 
equally reached by assisted adiabatic processes, such as 
those realized by STA. Therefore, reaching the minimal 
fluctuations in work does not necessarily require a work 
protocol to be slow. In addition to providing new insights 
into fluctuations phenomena in small systems, the results 
of this work should be of interest to the design of reliable 
and efficient energy devices at nano and micro scales. 


Appendix A: Ensemble-based quantum optimal 
control theory 

The quantum optimal control theory (OCT) we con¬ 
sidered in order to suppress the fluctuations in e~^ w is 
somewhat different from a traditional case due to two as¬ 
pects. First, we need to handle a thermal ensemble (see 
Eq. § instead of a single quantum state as the initial 
state. Second, the quantity we need to optimize is based 
on two-time energy measurements, which is not an ob¬ 
servable. Because of these peculiarities, it is necessary to 
outline some details in our OCT calculations. 

First, to suppress the fluctuations in e~@ w , an addi¬ 
tional control field is considered, with the total Hamilto¬ 
nian H of the system given by 

H(t) = H 0 (t) + H OCT (A(t)), (Al) 

where H OCT is the control Hamiltonian and A(t) is the 
time-dependent amplitude of a control field. The time- 
evolving state | (f>i) obeys the schordinger equation: 

I Mt)) = (A2) 

where | cf>i) denotes the i-th state in the initial thermal 
ensemble in the energy basis. 

Consider next a certain quantity described by Eq. ([3| 
as a target control function, denoted as L±. The job is 
to minimize the following: 

N 

pM\Mr))\ 2 nE l7 E'), (A3) 

hi=i 


where |0 ,(t) = U\tpi) is the final state evolved from the 
initial state | ipi) (eigenstate of H 0 ). For a control prob¬ 
lem, typically a cost function is also needed to reflect a 
cost-related constraint. This cost function can be con¬ 
structed as 


1,2 = \f «A 2 (t)df, (A4) 

where A(t ) is above-mentioned amplitude of the control 
field and k is a weightage factor. The overall target func¬ 
tion can then be defined as J = L\ + L 2 . That is, the 
problem is now to minimize J under the general dynam¬ 
ical constraint reflected by the Schrddinger equation. 

To proceed we introduce N Lagrange multiplier vectors 
as a function of t, denoted by {\h(t)),i £ N}. We then 
minimize J instead, with 

N T 

J = L\ + L2 + + i{h\H\(f>i(t))]dt. (A5) 

i=i J 0 

Let |(5be the variation in |</>;(£)) due to SA(t), an 
arbitrary variation in A(t), then the variation in J due 
to SA(t) is found to be 


N 


83= Y, PiMj\8Mr)) + h.c.]f(E i ,E , j ) 

i,j =1 
N 

+ J2[33 t )\ 8M t )} +h.c.} 

i =1 

N r 

[{k{t)\6<j>i(t)) + h.c.]dt 

i=l J 0 

N n. T 

+ i[{li(t)\H\6(j>i(t)) - h.c.]dt 

i=l J 0 

N rT arr 

+ X!/ 0 - h.c.]8Adt 


(A6) 


nA(t)SAdt. 


To minimize J we let 8J = 0. Since the variation is 
arbitrary, one has the following relations: 


' I U{t)) = -iH\k(t)) 

2Im[(k(t)\ — |&(t)>] + «A(i) = 0 

N 

p i J2W j )f(Ei,E' j ) + \l i (r))=0. 

. i=1 


(A7) 


The above list of relations can be numerically solved by 
an iteration procedure [29]. In our actual calculations, in 
precisely the same manner as Ref. 7], we introduced a 
time dependence to the “penalty factor” n to ensure that 
the OCT control field is zero at t = 0 and at t = t. 
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Appendix B: A classical version of the minimal work 
principle 

An extension of the minimal work principle proven in 
Ref- [23] to the classical domain can be carried out, fol¬ 
lowing essentially the same steps used in our above proof 
of the minimal classical work fluctuations. With the same 
notation as in the main text, the average work for a gen¬ 
eral process starting from Hq and ending with H is 


jR j H 

(W) + ( H c 0 ) = J [dI 0 } HS po(I 0 ) J dd 0 J [dI] HS H°(I) J dd 5[Ir - me T - 0\. 

Upon integration parts, we have 


jR 

(W) + ( H c 0 ) = J [dI 0 } HS p 0 (/o) J dd o J dd S[d T - d] 

| H c t (I R ) - f iL d [m(I)] H , o £ [dI') H c S[I T - I'] |. 


The parallel result for an adiabatic process is 


(W) + (H£) = J [ dI 0 ] H § PoVo) J dd 0 J dd 6[6 r - d] 


H C r(I R ) - £ d [HZ(I)] h§ £ (dI'] H§ S[I 0 - I') \. 


One hence finds the difference between (W) and (W): 


(Bl) 


(B2) 


(B3) 


(W) ~{W)= d[H c T (I)W 0 ©(/), 

Jl L 


where 


jR ( j 

0(/) =J [dI 0 \ H§ Po(Io) j ddA- j dd 6(dr - d) J [dl'] H§ S(I T - /') 


+ j dd S(d r - d) J [di'] H §6{io -/') >■ 


(B4) 


(B5) 


Note that due to the same-ordering condition, d[H!)(I)] H o is not negative. Splitting the integration in Eq. (B5) over 
[d/y] h c into two intervals: from I L to I and from I to I R , we have 


0(7) = / [dI 0 ]H§ 6 f dd 0 | [ dd S(d T - d) - / dd S(d T -d)[ (dl '] Hg 8(1 r - I') 


!l L 


r 1 e -PHS(io) r r r 1 

J [dioUs — Y* — J de ° J de s{8t ~ e) J IL [dI ' ]H o ~ J ') 


(B6) 
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For similar reasons as we derive Eq. (431, we now have the following inequality, 

e -PHZ(i) r i ( ■ 


0 (/)> 


z o Ji L 

e -0HZ(I) ,/ 


[dloUs dd 0 { / dd 6(9 T -B)~ dff 6(d T - d) / [d/V= <S(/ T - /') 


//*■ 


7 C / 
^0 JI 

e -f>H-(I) r-1 


Z 0 Jl L 

e-PHS(i) ,/ 

^0 7/r- 

e -0HZ(I) ri 


J 0 Jl L 


[dI 0 ]H 8 j ddo J dd 6(d r - d) f [dI'] H8 S(I T - 1 ') 
[<tfo]ffg j ddo J dd 5(9 T - 0) 

[d/oUg J dd 0 J dd5(d T -d) f[dI'] H o S(I T — I') 

r e -l3HSli) r [i 

[dIo\ H§ J ddo - J dd Jjdl'Us 


(B7) 


= 0. 


This finally leads us to the conclusion that (W) > (W). 
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